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Abstract 

We propose a nonperturbative formulation of chiral gauge theories. The method 
involves a 'pre-regulation' of the gauge fields, which may be implemented on a lattice, 
followed by a computation of the chiral fermion determinant in the form of a functional 
integral which is regularized in the continuum. Our result for the chiral determinant is 
expressed in terms of the vector dike Dirac operator and hence can be realized in lattice 
simulations. We investigate the local and global anomalies within our regularization 
scheme. We also compare our result for the chiral determinant to previous exact £- 
function results. Finally, we use a symmetry property of the chiral determinant to 
show that the partition function for a chiral gauge theory is real. 
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1 Introduction 



The non-perturbative formulation of chiral gauge theories is a long-standing problem of 
quantum field theory with both practical and theoretical implications Q. Recently there 
has been a rekindling of interest, largely due to 'tHooft [0, in the idea of using gauge field 
interpolation to couple lattice gauge fields to continuum fermions. By keeping the fermions 
in the continuum one avoids the difficulties associated with realizing chiral fermions on the 
lattice ||. This idea was discussed previously in the literature in j|, ||, [| and has now 
recently been further developed in |7], ||, |jj . 

In this approach, the gauge fields and fermions are treated differently, with the gauge 
fields originally defined only on the links of a spacetime lattice, as in the usual formulation of 
lattice gauge theory. However, an interpolation algorithm associates the gauge link variables 
U^(x) with a continuum gauge field A^x) in whose background the fermion part of the 
functional integral is evaluated. Thus, given an interpolation scheme the problem that 
remains is to give a continuum formulation of the chiral determinant which can be evaluated 
to a desired accuracy within a finite computation. 

In this paper we give a simple, gauge invariant formulation of the chiral determinant in 
terms of the vector-like Dirac operator which can be realized, at least in principle, in lattice 
simulations. The relevant background fields for our determinant have been 'pre-regulated' by 
the lattice interpolation and for the purposes of our analysis we will take them to be smooth, 
with variation on length scales larger than a chosen scale a, where a is related to the lattice 
spacing a. We will discuss below to what extent various interpolation schemes satisfy this 
property. At the end of the procedure we allow a, a — > 0, but only after first taking to 
infinity the continuum mode cutoff iV used in the functional integral. The good behavior of 
the background field allows us to make well-defined manipulations of the functional integral 
and in particular to separate low frequency physics from the high frequency physics which 
comes from modes near N. 

Our partition function is 



where the sum is over all gauge link configurations, the action Sym is the usual Yang-Mills 
lattice action, and the determinant is a functional of the continuum gauge field A^x) which 
is uniquely determined from each discrete set of links {U^(x)}. The determinant det[ JpL } of 
the chiral Euclidean Dirac operator in the background field A^(x) will be defined in section 
3. As we will discuss below, it is also straightforward to define regulated fermion correlators 
within our scheme. 

Let us recall some well-known results concerning fermion determinants. In Euclidean 



Z = e~ s ™W det[p L ], 




{u} 
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space fermion determinants for vector-like models are real and positive semi-definite. It is 
possible to convert a chiral model with fermions in representations r into a vector-like one 
by the addition of 'mirror' fermions which are exact copies of the originals, but in complex 
conjugate gauge representations r*: 

det[ pv } = det[ p L r } det[ p L r * ] = det[ fl L r ] det[ p L r ]* . (1.2) 

Thus, we see that the magnitude of the chiral determinant is simply equal to the square 
root of the corresponding vector-like determinant. We also see that any gauge anomaly in 
the chiral determinant must be a pure phase, since a vector-like model is anomaly free. In 
this paper we are primarily interested in models in which gauge anomalies are absent, so the 
important quantity is the non-anomalous phase of det[ Pl ], which reflects the chiral nature 
of the model. 

Our regularized chiral determinant, described in detail in section 3, is given in terms of a 
regularized fermion functional integral. We regulate the integral by truncating the functional 
measure to a finite number of modes, N. This cutoff renders the integral finite, but at the 
cost of introducing some gauge non-invariance in the magnitude of the determinant. Since, 
as mentioned above, we already know the desired magnitude of the fermion determinant in 
terms of the corresponding vector-like model, this does not lead to any ambiguity in defining 
the regularized model. We simply modify the magnitude to agree with the square root of 
the corresponding vector-like model. 

We will argue that our determinant has the following properties: 

(1) It correctly reproduces the known global and local anomalies, as well as allowing for 
anomalous phenomena such as fermion number violation. (This is in some sense 'built-in' 
to our formulation as it mimics the naive continuum definition as much as possible.) 

(2) It yields a gauge invariant result in the limit N — > oo, at least for infinitesimal gauge 
t r ansf or mat ions . 

(3) Its construction is accomplished using only the eigenvalues and eigenfunctions of the 
vector-like Dirac operator p. 

This paper is organized as follows. In the following section we briefly discuss some 
interpolation schemes and their properties. In section 3 we give our definition of the chiral 
determinant in terms of a regularized fermion functional integral. In section 4 we investigate 
the local and global anomalies within our regularization scheme. We show that in models in 
which gauge anomalies cancel the phase r) of our determinant is invariant under infinitesimal 
gauge transformations even for finite truncation to N modes. In section 5 we discuss how 
our results are related to previously derived exact representations of the imaginary part 
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of ln(det[ Pl })■ In section 6 we discuss a simplification of the partition function which 
arises from the behavior of the chiral determinant under reflection of the background field. 
We conclude with a summary and an appendix on the convergence properties of the chiral 
determinant. 



2 Pre-Regulation of Gauge Fields 

In this section we give a brief overview of some possible interpolation schemes, with em- 
phasis on their smoothness properties. As we will see in later sections the required size of 
the continuum cutoff on the fermion modes depends on the smoothness properties of the 
background gauge fields. It is desirable that the interpolated fields have their support in 
momentum space concentrated at momenta less than some scale 1/a, which is presumably 
controlled by and of order the lattice spacing a. It is also desirable that the interpolation 
scheme be gauge covariant, so that the effect of a lattice gauge transform on the links is 
consistent with effect of the interpolated gauge transform on the continuum field. More 
explicitly, under a gauge transform we should have 

{U} -> {U*} 
A,(x) Af{x), (2.3) 

where Q is the lattice gauge transform (valued only on discrete lattice sites) and Q c its con- 
tinuum interpolation. Even given the above criteria the choice of interpolation prescription 
is highly arbitrary. 

The interpolation scheme based on the geometrical definition of topological charge [11 



obeys the gauge covariance condition, but both it and the scheme based on minimizing the 
Euclidean action proposed in produce continuum fields which are in general only piece- 
wise continuous. Without some additional smoothing, these configurations possess Fourier 
transforms with support that falls off at large momentum k only like one over \k\ to some 
power. In general we would prefer configurations which are completely smooth (infinitely 
differentiable) and hence by the Riemann-Lebesgue lemma have Fourier transforms that fall 
off exponentially at large momenta. 

A very simple scheme which leads to smooth continuum fields was recently given by 
Montvay ||. (Unfortunately this scheme is specific to U(l) gauge theories.) First one 
imposes a gauge fixing such as Landau gauge to define a set of gauge fields A XfJ/ on the 
discrete spacetime lattice. (This can always be done in such a way that the A Xfl are bounded 
by 2ir/a.) The Fourier coefficients A^ of the lattice gauge field A xtl are valued only on 
discrete values of k^ in the Brillouin zone: k^ = 2Tm^/L^ , < < L^/a. (L^ is the 
extent of the lattice in the /i direction.) One can use these Fourier coeffi cents to define a 
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continuum gauge field through 



M*) = irxfe E^'^V (2-4) 

The resulting continuum gauge field is infinitely differentiable, and has zero support in 
momentum space for \k\ > 2ir/a. A slightly modified version of this scheme can be given 
H which also satisfies the covariance condition ( |2.3|) . An interpolation like (|2.4j ) relates the 
lattice and continuum gauge transforms and guarantees that the latter also have support 
only in the Brillouin zone. 

For our purposes in what follows, we will assume that the gauge field backgrounds in 
which our functional determinant is to be evaluated are bounded and infinitely differentiable, 
and hence that their Fourier transform falls off exponentially rapidly at momenta large 
compared to a scale 1/a which is controlled by the lattice spacing a. We will also restrict 
ourselves to smooth gauge transforms Q(x), such as would result from an interpolation of the 
type above from a lattice gauge transform. This implies that in the band-diagonal matrices 
which will appear in the next section, the entries which are outside the bands can be made 
zero by choice of the width of the band. In the specific case of an interpolation like that of 
|J, the width of the bands will be ~ 1/a with entries outside the bands exactly zero. 

We should note that while this assumption of smoothness of the background field allows 
us to discern various nice features of the determinant, such as the band-diagonal properties 
of certain matrices, it is not absolutely necessary to demonstrate the three main properties 



listed in the introduction. For that purpose the interpolations of |1CJ, [11], |J with their 
piecewise continuity should suffice. 



3 Chiral Determinant 

Our task is now to define the functional determinant in the smooth background A^x). We 



will follow a straightforward procedure similar to that first used by Fujikawa ]T2[ in his func- 
tional integral approach to the anomaly. We work in Euclidean space-time which results from 
the Wick rotation x° — > — ix 4 and Aq — > iA^. The Dirac operator If) = ^^D^ = 7 M (<9 M + y4 At ) 
becomes a Hermitian operator 

p = 7 4 D 4 + Y D ., (3.1) 

We use the convention that 7 is Hermitian and 7 fc anti-Hermitian. The Hermitian 7 5 is 
defined as 

75 = = -7 1 7 2 7 3 7 4 - (3.2) 

The Euclidean metric is g^ u = (—1, —1, —1, —1). 
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The determinant is defined, formally, through the fermion functional integral 

det{ p L \ = J Dtp i Dtfj L e'/^ lWi . (3.3) 

The chiral fields are given in terms of the complete set of orthonormal eigenfunctions of the 
Hermitian Dirac operator: 

p<P n = A n n . (3.4) 

As usual, we imagine that our system has been placed in a box of size L with appropriate 
boundary conditions in order to discretize the eigenvalues. We define a complete set of chiral 
modes by 

<P L n = V2P L( j) n A n >0 (3.5) 
= Pl4>u A n = 0, 

where Pl = {—^-)- The modes 0^ are defined equivalently but with left-handed projectors 
replaced with right-handed ones. The basis {0^, 0^} is complete and orthonormal, as we 
can see since pj 5 (p n = — A n 750 n - 

A subtle but important point should be emphasized here: the choice of positive eigenvalue 
modes which are used to define {0^, 0^} must be determined for some fiducial value of A^(x), 
which we will take to be A^(x) = 0. As the gauge field is varied from zero to some arbitrary 
configuration A^(x), the evolution of modes which initially had positive eigenvalues A n > 
must be followed in order to define the new set of modes in the A^x) background. In the 
process, some of the initially positive eigenvalue solutions may end with negative eigenvalues. 
It is this phenomena that is responsible for possible sign changes of the chiral determinant^. 

If we had not followed the procedure of 'tracking' the modes from the fiducial background 
to the background of interest, but rather applied the definition (|3.5|) in a naive way at each 
value of A^(x), discontinuities in functional derivatives such as SA S ^ det[ pL ] could arise. 
This sign ambiguity in the definition of the determinant is the cause of the well-known global 
anomaly, as first described by Witten |0| 0. We will return to this point in the next section 
when we show that our regularization of the determinant reproduces the correct results for 
the global anomaly. 

Having chosen an orthonormal basis for the background, we can then expand 

MX) = E bn<Pn(xV- (3-6) 

Hn practice, this tracking of eigenmodes would be a rather cumbersome procedure, requiring an inter- 
polation of the background field of interest to (x) — and the computation of the low- lying eigenmodes 
over this interpolation. 

^Fujikawa fll^ is primarily interested in the effect of infinitesimal transformations on the functional 
measure. His formulation, which ignores the tracking of eigenmodes, is sensitive only to local anomalies. 
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With this expansion, the functional integral takes on a particularly simple form 
det[ p L ] = / II II db n da m e i S„ A ^ a ™ det{ C L ] ■ det[ C R } 

n m 

= det[ C L } ■ det[ C R } J] i \ n . (3.7) 

n 

Fermion correlators can also be explicity evaluated, for example, the propagator 
(Mx)MO)) = J Dvb L Dij> L (Mx)M0)) e~I d4x ^^ 

fe t 0^)^(0)) det[ C L } ■ det[ C R } ft i X n . (3.8) 

\ m / n 

The factor Y\ n X n is simply the square root of the corresponding vector-like functional integral 
(i.e. with no chiral projector in the action of ( |3.3| ) ), up to the sign ambiguity we discussed 
previously. The factors of i in the infinite product amount to an overall phase factor which 
is A^(x) independent and are irrelevant to our analysis. We will drop them in what follows. 

The additional Jacobian factors det[ C L } , det[ C R } arise from the change in basis we 
have had to make from an initial fiducial basis which we define as a product over free modes. 
In other words, we choose an initial fiducial measure for the functional integral 

db n da m , (3.9) 

n m 

where the coefficients b n , a n are those that arise in the expansion of ( |3.6| ) in terms of free 
solutions. In order to perform the integral using the orthonormality properties of solutions 
in the A^(x) background, we have to change basis and hence the extra Jacobian factors must 
appear. 

The matrices C L , C R are defined as follows. Let the lack of an additional superscript 
denote free eigenmodes and the superscript A denote eigenmodes in the A^x) background. 
Then 

= J d*x <^(*M(z) 

= J d 4 x <Pi(xY P L n . (3.10) 

Similarly, 

cL = (0*10^) 

= J d*x <P R (xY^ A (x) 

= J d 4 x <f> n ( x y p R <j>i{x). (3.ii) 
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The C matrices are complex but unitary, so the factor det[ C L ] ■ det[ C R } is formally a 
pure phase. This phase, when combined with the potential sign changes in Y\ n X n (equivalent 
to phases of wr), constitute the chiral phase information mentioned in the introduction. In 
the absence of gauge anomalies, det[ pi } should be gauge invariant, which in turn requires 
that the result (|3.7|) is gauge invariant. At the formal level, a gauge transformation acts as 
a unitary transformation on the C L,R matrices, and hence should leave det[ C L ] ■ det[ C R ] 
invariant. We will examine the gauge invariance properties of our regulated version of this 
object in section 4. 

In what follows we will regularize all of our previous expressions by eliminating all but a 
finite number of eigenmodes from our functional integral. This will truncate all infinite sums 
and products to finite ones, and also infinite matrices to finite matrices. All expressions and 
manipulations will then be well-defined and finite, and could in princple be implemented on 
a gedanken-computer. (See |7|, |i| for work on similar regularization schemes.) Of course, 
this regularization is precisely a 'hard cutoff' in the space of eigenf unctions, and violates 
gauge invariance. One of our main results will be that for pre-regulated backgrounds A^x), 
the violations of gauge invariance are limited and readily compensated. In particular, gauge 
non-invariance due to the finite truncation will be seen to only affect the magnitude of the 
chiral determinant, allowing the chiral phase information to be extracted in a well-defined 
manner. Our truncated expression for det[ pi ] is then 

N 

det[ C L ] ■ det[ C R } J] X n , (3.12) 

n 

and for the propagator 

/ N \ N 

(M*)M0)) = Ei«(0) det[ C L ] ■ det[ C R } J[X n . (3.13) 

\ m / n 

where the C matrices are now finite dimensional. The number of modes that are kept is 
N ~ (LA) 4 , where A is roughly the UV scale associated with our mode truncation. (This 
is up to additional factors due to degeneracies and zero modes.) Here L is the size of our 
box. As we will specify, the scale at which our truncation is made is determined by the 
smoothness scale a, the lattice spacing a and the box size L. 

The pre-regulation (smoothness) of the gauge field implies that the C matrices have a 
very simple form. We can see this by first making a useful observation about the explicit 
forms of the modes which appear in (|3.10|) and ( 3.11|) . Because the free basis is essentially a 



plane wave basis, it is useful to examine our eigenvalue equation (p.4|) in momentum space. 
The eigenvalue equation for 4> n (q), has the form 

{ii-X n )4> n {q) + / 'd 4 k 4,(q-k)<p n (k) = 0. (3.14) 
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Here A^k) represents the Fourier transform of the background field A^(x). Because of the 
smoothness property of A^x), the Riemann-Lebesgue lemma tells us that A(k) — > as 
\k\ —>■ oo. In fact, A(k) goes to zero exponentially rapidly for \k\a » 1. 

We will now show that the solution to (ft.!4|) , <j) n (q), has its support only in regions of 
momentum space centered around values of q which satisfy q 2 = \\. The size of those 
regions of support is of course determined by the properties of A^(x). To simplify ( p,14|) , let 



us choose the Dirac basis for our gamma matrices so that 74 is diagonal, and the rest are 
off-diagonal. Now let us choose a frame in which the momentum g« = (0, 0, 0, q^). We will 
show that if q^ is sufficiently different from ±A n , there is no solution to (gig) . 

Let the the four-spinor (j) n (q) have the two-spinor components u n (q),v n (q). In the frame 
we have chosen, the eigenvalue equation becomes 



(<?4 + K)un(q) = I d 4 k (iA A (q - k)u n (k) + A { (q - k)uiV n {k)) (3.15) 
(<?4 - K)vn(q) = I d 4 k (iA 4 (q - k)v n (k) + A^q - k)<JiU n (k)) . (3.16) 



Multiply the top equation in ( |3.15| ) by u^q) and the bottom by v^q) and integrate both 
over a ball B of size a -1 centered about our chosen q* = (0, 0, 0, qX) (so q = q* + q', and 
we integrate f B d 4 q' with \q'\ < a -1 ). 

We can then rewrite the equations as 

(qt - A n ) = ••• 

(<Z 4 * + An) = • ■ • , (3.17) 

where the terms on the right hand side denoted by • ■ • are bounded, as we will see below. 
This implies that there is no solution when \q\ ± A n | are both taken sufficiently large. 

The terms on the right hand side are of three types. In matrix notation, the first type is 
of the form (we suppress numerical factors like 7r) 



(u 


A 


u) 


(u 


u) 



< (L/a) 2 a-\ (3.18) 

To compute the bound in ( |3.18| ) we have used Parseval's theorem, 

Jd A k\A(k)\ 2 = Jd 4 x\A(x)\ 2 < L 4 a- 2 (3.19) 

and the result that A(k) has support only for momenta \k\a ^ 1. The second type of term 
is 

\u\u) 

where the bound applies if u n and v n both satisfy ( |3.15| ) and for q\ sufficiently different from 



both ±A n . (If u n and v n were left arbitrary, one could choose \u n \ — > while v n remains 
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fixed, which would make the left hand side of ( |3.20| ) arbitrarily large. However, it is easy to 
see that such u n and v n cannot satisfy ( [3.15 ) when the lhs of those equations are sufficiently 
large.) 

The last type of term is of the form 



J B d 4 q' q' ul(q)v n (q) 



(3.21) 



Js^Y ul(q)u n (q), 

and is bounded by a -1 if again both u n and v n satisfy ( |3.15|) and q\ is sufficiently different 
from both ±A n . 

Thus we conclude that eigenfunctions which satisfy ( |3.14| ) have support only in regions 
of momentum space which lie within a distance of order 



max[ (L/a) a 



^ 2 -- 1 a- 1 



(3.22) 



of values q which satisfy q 2 = A 2 . 

Intuitively, we can understand this result as follows: the smooth background A has 
support in a compact region of momentum space, and is bounded in magnitude as well 
(\A\ < a^ 1 ). Its effect on eigenmodes is to 'mix-up' the original plane wave modes, but only 
those that are within a certain band of each other, whose size is determined by A. The 
overlap between modes <f) m , <p n is zero for sufficiently large \m — n\. We define the number 
Na such that the overlap between any modes with \m — n\ > Na is zero. 

The C matrices therefore have the band-diagonal form displayed below. We will always 
assume that the size of the matrix N is much greater than the width of the band Na- This 
roughly corresponds to a choice of cutoff A for the regularization of our determinant. We 
note that from ( |3.22| ) it is clear that A must always be much larger than the corresponding 
lattice cutoff a -1 . 
























Vo 

































o\ 
























(3.23) 
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For simplicity, we have adopted an index in the above matrix representing 'momentum' 
k rather than the index n. The latter runs from to oo whereas the momentum can be 
positive or negative. We are pretending that there is only one component of momentum - in 
reality the C matrices are actually multidimensional with additional labels representing the 
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individual momenta ki, as well as internal group indices. The multidimensional matrices have 
their band structure centered about vectors ki in momentum space satisfying J2t=i kf ~ 
for some A n . 

Now we can see the problem that arises with truncation to a finite number of modes: the 
finite dimensional matrices C L , C R are no longer unitary. In fact, the products C L C L ^ and 
C R C R * are no longer the identity but have the following form: 



\ 



V 



(3.24) 
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Let us denote the non-diagonal submatrices in the upper left and lower right generically as 
X L,R . They are of dimension Na and are the result of a loss of unitarity that comes from 
our truncation. Since an infinite number of eigenfunctions are necessary to span the space 
of solutions, unitarity requires an infinite sum: 



fimp — E C mn C n p 



oo 

E 



t A \ 



) 



t\rt A ) 



(3.25) 



From previous analysis we know that overlaps between modes m and n for which \m— n\ > Na 
are negligible. Therefore, for m,p < (N — Na) the unitarity relation ( |3.25|) is unaffected. On 
the other hand, for (N — Na) < m,p < N, intermediate modes in the sum with nontrivial 
overlap are removed in the truncation, and hence the X L,R matrices in the figure above are 
no longer necessarily close to the identity. Note that the dimensionality of this matrix is 
independent of the parameter N and the mode truncation in general as long as iV >> iV^. 
From the form of (|3.24|) , we see that the product det[ C L } ■ det[ C L ^ } is now equal to: 



det[ C L } ■ det[ C if ] = det[X L ], 



(3.26) 



and similarly for the C R matrices. The determinants of X L,R are real, as X L,R is Hermitian: 

N 



(■^ )mp 



n=N —N a 
N 

n=N-N A 

(XI) pm = {X ) rapi 



(3.27) 



10 



with similar result for X R . 

As noted previously, the formal unitarity of the infinite dimensional C matrices implies 
that det[ C L ] • det[ C R ] is a pure phase. For finite dimensional C matrices the constraint 
that remains from unitarity is that 

| det[ C L } ■ det{ C R } | 2 = det{ X L } ■ det{ X R ] = det[ X ] 2 , (3.28) 

which allows for an arbitrary phase in det[ pi ], but does not restrict the magnitude to be 
unity. Here we define X as a diagonal matrix diag{Af, • ■ • , A^ A } , where the \* are given 
by the square root of the product of the corresponding n-th eigenvalues of X L ,X R . X is 
Hermitian and det[ X } is real, so we can still extract the phase unambiguously from the 
finite dimensional C matrices, where 

det[ C L } ■ det[ C R } = e ir > det[ X } (3.29) 

or 

rj = kn (in (det[ C L ] ■ det[ C R ])) . (3.30) 

In the next section we will examine the behavior of det[ C L } ■ det[ C R ] under gauge 
transformations. Any anomaly, being a pure phase, resides in the factor e ir> . In the absence 
of gauge anomalies, r] is gauge invariant, which is equivalent to the requirement that under a 
gauge transformation, the change in det[ C L ] • det[ C R ] is purely real. We will see explicitly 
that this is the case in section 4. 

Our final result for the regularized chiral determinant is 

N 

det[ Jp L } reg = e*" J] X n , (3-31) 

n 

where the phase rj is gauge invariant in the absence of gauge anomalies and must be extracted 
from the finite dimensional C matrices via ( ET29D or flOOl) . In the append ix we discuss the 



convergence properties of n. It appears that rj is determined by the Na x A^^ submatrices of 
C L,R (i.e. it depends on the properties of eigenmodes with n ^ Na) , and therefore converges 
to a well-defined value as — *• oo. An alternate method of extracting the complex phase in 
( p.31| ) already exists, as we will discuss in section 5. 



4 Local and Global Anomalies 

In this section we examine the local and global anomalies within our regularization scheme. 
We will verify that our regularization scheme leads to the usual anomalous Ward-Takahashi 
(WT) identities in the limit that A^ — ► oo and that in the absence of anomalies the phase fac- 
tor rj is gauge invariant. Finally, we shall examine how Witten's global anomaly is manifested 
in our scheme. 
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First let us review the behavior of the functional measure under rotations of the fermion 
fields All: 



i> L -> $ L e ia{x) . (4.1) 

These rotations can also correspond to non-Abelian gauge transformations if we allow a(x) to 
be an element of the Lie group: a(x) = a a (x)T a . Any subtle effects from such a transforma- 
tion are to be found in the functional measure, or equivalently in the Jacobian determinants 
det[ C L } ■ det[ C R }. As usual, for infinitesimal rotations we can write 



det[ C L ] = det[ 5 mn + i J d A x a(x)(p^ n (x) J( <f)^(x)} 

= exp[ i J2 I d A x a(x)ri(x)W(x)] (4.2) 

n J 

and similarly, 

det[C R ] = exp[-iJ2 f d 4 x a{x)<f) R {xy<j) R {x)]. (4.3) 
Combining these equations gives 

det[ C L } det[ C R } = exp[tJ2 f d 4 x a{x){^{x)^ (j)^{x) - <j) R {x) ] (j) R (x)) } 

= exp[i^2 J d 4 x a(x)<i) n (x) ] ^(t) n {x)}. (4.4) 



In Fujikawa's scheme | 12fl , the infinite sums are regulated by the insertion of a convergence 
factor of 

/(^/A 2 ) = f(X 2 JA 2 ), (4.5) 

where /(0) = 1 and /(oo) = /'(oo) = /"(oo) • • • =0. The result is unchanged as long 
as the function / is smooth and obeys the above boundary conditions. The usual choice for 
/(^/A 2 ) is /(^/A 2 ) = exp(-p>/A 2 ). 

Using Fujikawa's result in the limit N — > oo (which corresponds to taking A — > oo) yields 
the well-known result 

oo 

lim £0„(*)t 75 0„(x)/(A 2 /A 2 ) = -^FF(x). (4.6) 

A— >oo u 
n 

If the original transformation ( |4.1| ) had been a non-Abelian gauge transformation, (|4.6| ) 
would have additional color structure and be proportional to Tr[{T a ,T b }T c }. 

One could instead have taken the function f(]ft/A 2 ) to approach a step-function, so as 
to reproduce a hard mode cutofff (jTj], |7j . Since the result is independent of the detailed form 
of /, we obtain the usual anomalous WT identity for the regulated current 

<V M = ^ FF(x). (4.7) 
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Note that the fermionic currents J^{x) are themselves divergent objects (they involve a 
product of operators evaluated at the same spacetime point x) which require regularization 
and some choice of subtraction in their definition. In our scheme any correlator is regulated 
automatically and non- locally by the mode truncation. For example, see ( |3.13| ). As x — > 
the usual short-distance divergence is cut-off by the truncation of the series at N. 

From ( f4.4j) we see that the absence of anomalies implies that the phase factor rj defined 
by the N — > oo limit of ( |3.30|) is gauge invariant, at least under infinitesimal rotations. The 
corrections that would appear on the RHS of ( |4.6|) due to a finite truncation are of the form, 
e.g., TrF k /M 2k ~~ 4 (k > 2) and vanish as N is taken to infinity. What is more, the corrections 
correspond to local operators and can be compensated by proper choice of local counterterms 
if desired. 

Let us now consider the global anomaly flT3"H , which arises in SU(2) gauge theories due 
to the fact that the fourth homotopy group of 577(2) is nontrivial, 

n\SU(2)) = Z 2 . (4.8) 

This means that in four dimensional Euclidean space there is a gauge transformation Q(x) 
such that Q(x) — * 1 as \x\ — * oo, and Q(x) covers the gauge group in such a way that it 
cannot be continuously deformed to the identity. The so-called mod two Atiyah-Singer index 
theorem then implies that as a gauge background A^(x) is changed continuously to its 
gauge transform A^(x) (e.g. via A l (x) = A^(x) + (1 — t) (A^x) — A^(x)) as t goes from 
zero to one), an odd number of positive- negative pairs of eigenvalues of pwill switch places. 
(The set of eigenvalues must match exactly at t = and t = 1 since A^(x) and A^(x) are 
related by a gauge transform.) This leads to an overall change in sign in the product 

oo 

II X n (4-9) 

n 

as long as we choose our eigenvalues in the way we have described in the previous section, 
tracking modes continuously beginning from some fiducial background A^(x). 

Now, given that we are only interested in backgrounds which arise from our pre-regulation 
procedure, it is easy to see that any modes which switch places are within roughly Na of 
n — 0. Therefore, as long as >> Na, the truncated product 

N 

II X n (4-10) 

n 

which appears in our regulated determinant undergoes the same sign change as the infinite 
product above when the background field is transformed by Q(x). (Of course, this Q(x) 
must relate two backgrounds A and A' which satisfy our smoothness conditions, so Q(x) 
itself must be smooth.) The modes near the cutoff, — Na < n < N, are not necessary 
to reproduce this result. 
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5 Exact Representations 



In this section we compare our form of the determinant to exact representations previously 
obtained using ^-function methods || 15| . We will not give the derivation of the results of 



|5|, [L5[] here, since the details are somewhat technical, but will merely state them. Suppose 
we wish to compute the imaginary part of the chiral effective action ln(det[ pi ]) in the 
background A^x) relative to ln(det[ pi ]) in the fiducial background A^x) = (we take 
the latter to have no phase as a choice of convention). First form the five dimensional 
background gauge field A t which interpolates adiabatically between A^x) = and our 
chosen A^x) as the parameter t varies from — oo to +oo. Next, consider the 5-dimensional 
Dirac operator defined by 

$> = ( *l 5 d 5 + pi ), (5.1) 

where pi is the vector-like Dirac operator in four dimensions. The result in an anomaly free 
model is the following: 

Im ln{det[ p L ]) = n ( 77(0) + dim ker Jfc ), (5.2) 

where 7](0) is the famous '^-invariant' of the Atiyah-Singer Index theorem. It is given by the 
analytic continuation to s = of 

V(s) = £i^> (5.3) 

A^O 

where A denotes eigenvalues of p$. One can think of 77(0) as the regularized spectral asym- 
metry of p^\ 

V(0) ~ El" E 1 - (5-4) 

A>0 A<0 

The term dim ker p§ simply counts the number of zero modes of the operator p§. We 
see that the determinant changes sign whenever such a zero mode appears. An identical 
change in sign can be seen to result in our treatment if we recall Witten's result [13], using 



a construction like the five dimensional one above, that a zero mode of p$ implies that in 
the interpolated background A t the flow of eigenvalues of p^ is such that an odd number 
of eigenvalue pairs change sign. Thus the sign change in due to dim ker in the exact 
representation matches that due to eigenvalue flow in our form of the determinant. 

This leaves the phase rj defined in ( |3.30|) to be identified with 77(0) . This identification 



is extremely nontrivial mathematically, as it relates the ^-invariant in five dimensions to 
some rather detailed properties of the eigenfunctions of the four dimensional Dirac operator. 
As mentioned previously, we have not rigorously proved (although it is plausible - see the 
appendix) that our phase rj converges to a well-defined limit as the cutoff is taken to infinity. 
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On the other hand, r](0) is defined in terms of an analytic continuation which, while math- 
ematically sound, may or may not have the same physical content as our 77 0. These issues 
clearly deserve further investigation. 

6 Reflections and Phases 

In this section we make an observation which simplifies the form of the partition function 
( |1 . 1|) . We show that despite the possibility of complex phases in the chiral determinant, in 
the absence of gauge anomalies the partition function itself is real. Our main observation is 
that the 77(0) part of the fermion effective action changes sign when the gauge background 
is reflected through any of the hyperplanes: {x M = 0} - in other words, it has odd parity. 
This implies that all imaginary parts eventually cancel in ( |1.1Q , leading to a real partition 
function. 

Consider a background A^(x) and its reflection through x 4 = 0: 

A*(xi,x 4 ) = Apfa, -x 4 ). (6.1) 

The five dimensional interpolations t A fl (x) and 'A* (x) will also be reflections of each other. 
We will specialize to the temporal gauge: A4 = 0. We do not lose any generality by doing so 
if we are working in a model without gauge anomalies, since in that case det[ pi ] is gauge 
invariant. Note that the pure gauge actions SVm(/4] and SVm|/4*] are identical. 
We want to show two things: 

(I) The phases of the determinant det[ pi } in the backgrounds A^(x) and A*(x) are related 
by a minus sign. In other words, 

7^(0)1 t A = - 7^(0)1 t A . , (6.2) 

which is equivalent to showing that there is a mapping of A n — > — A n when the gauge 
background is reflected through x 4 = 0. 

(II) The eigenvalues of p± are invariant under A — > A*. This should be clear since the 
Euclidean Dirac operator can be regarded as a Hamiltonian, and the energy eigenvalues are 
invariant under reflection of the gauge background. We will also see this explicitly below by 
looking at eigenf unctions and their eigenvalues. 

§In the work of Ball and Osborn ||T5[ similar results are derived using Pauli-Villars rather than (^-function 
regularization. 
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Given points (I) and (II), we can arrange the partition function in the following way: 

Z = J2 e~ SYM[A] det[ p L } 

{A , A*} 

= ]T e ~ SYM[A] 2 Re( det[ p L ] ) 

{A} 

£ e" s ™ [A] Re( det[ p L ] ), (6.3) 

{A , A*} 

where first and third sums are over all configurations and the second sum is only over the half 
of the possible configurations which remain after modding out by the Z2 reflection symmetry. 
We have also used the fact that the sign of det[ pL ] is the same in the A and A* backgrounds. 
This follows from (II) and the tracking of eigenvalues. 

Now to the proof of (I). The eigenvalue equation p^(p n = \ n (p n is as follows (we use the 
chiral basis for gamma matrices): 

[+ id 5 u n - id A v n + DiCTiVn] = \ n u n (6.4) 
[- id 5 v n - id 4 u n - Di(JiU n } = X n v n . (6.5) 

One can check that, corresponding to the original eigenfunction of p^ in the 1 A background, 

( U fX (6.6) 
is an eigenfunction of p^ in the background *A*, with eigenvalue — A n : 

fv n (Xi,-X 4 ,X 5 )\ 
\U n (Xi, -x A ,x 5 )J ' 

This shows that the eigenvalue spectrum of p$ in the reflected background is the negative 
of the original spectrum, which is sufficient to prove (I). 

We can see that this result also follows from the functional integral form of the deter- 
minant if we continue to work in the temporal gauge. In this gauge a parity-like operation 
relates the bases {<fin( x )} an d {0^(^)1 in backgrounds which are related by a reflection 
through one of the hyperplanes {x M = 0}. 

Given an eigenfunction (in the chiral basis where 75 = diag{l, —1}) of the Dirac operator 
in the background A: 



(*) = ( 1)1 I > ( 6 - 8 ) 



the following spinor function is an eigenfunction with the same eigenvalue, but of the Dirac 
operator in the background A*: 
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(Note that this proves point (II) above.) This similarity between the left and right chiral 
bases in the two backgrounds yields the following matrix relations: 

C L (A) = {C R (A*)}* 

C R {A) = {C L {A*)Y . (6.10) 

Together, these imply that 

det[ C L (A) } ■ det[ C R {A) ] = {det[ C L {A*) } ■ det[ C R {A*) ]}*, (6.11) 

which is equivalent to (I). 

We should note a limitation of the result: the insertion of an operator into the sum in 
( |6. 3D (e.g. to compute an n-point correlator) will in general destroy the reflection symmetry 
of the terms in the sum. Therefore we can only compute correlators using the weighting 
Re( det[ pi } ) if the operators are themselves invariant under at least one reflection, such 
as 

(0(xi,x 4 ) + 0(xi, -x 4 ) ), (6.12) 

where O is a generic field operator. This does not seem to be a significant limitation on 
our ability to investigate issues of interest such as chiral symmetry breaking, confinement or 
particle spectra. 

7 Conclusions 

We have proposed a formulation of chiral gauge theory which should allow, in principle, the 
evaluation of any quantity to a specified accuracy within a finite computation. In order to 
avoid the problems associated with lattice chiral fermions |§ , we employ an interpolation of 
the original lattice gauge fields to the continuum p], f|, [J. All of the fermionic aspects of 
the theory have been expressed here in terms of the eigenvalues and eigenf unctions of the 
vector-like Dirac operator, which can itself be implemented directly on a lattice if desired. 

The main focus of our investigation was the chiral determinant defined in terms of a 
regulated continuum functional integral. When the gauge field background is suitably well- 
behaved (e.g. resulting from lattice interpolations like those described in section 2), we find 
that an object with all the desired properties can be extracted from the functional integral 
as long as the mode cutoff N used to regulate the integral is kept sufficiently large. For well- 
behaved background gauge fields, any violations of gauge invariance introduced by a hard 
mode cutoff are confined to modes close to the cutoff. These violations of gauge invariance 
do not affect the complex phase information which characterizes the chiral nature of the 
theory. The order of limits that are necessary for our regulator are as follows: first, the 
continuum cutoff iV — > oo, followed by the original gauge lattice spacing a — > 0. 
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Our result for the chiral determinant is (from ( |3.31| )) 



det[ fl L ] reg = e*> J] *n , (7.1) 



where the phase rj is defined in ( |3.30| ) and the eigenvalues A n are described below (|3[ 
Both 7] (in the absence of anomalies and in the limit iV — > oo) and the eigenvalues are gauge 
invariant functions of the background field. In the appendix we give some evidence that 
suggests that det[ p L ] as defined above will converge to a well-defined value as N — > oo, 
with its value mainly dependent on modes with n ^ Na- ( |7.1|) is consistent with the exact 
representation derived using ^-function regularization, given the identification of the phases 
rj and 77(0) described in section 5. One difference between the functional integral formulation 
and the ^-function regularization is that in the former the physical aspects of the chiral 
determinant such as level crossing and the role of high versus low frequency modes are more 
transparent. 

Using our result for the chiral determinant we showed that despite the sum over complex 
phases, the partition function for a chiral gauge theory is real. The simplified result for 
the partition function that we derived in section 6 still requires the computation of the 
non- anomalous phase rj, as 

Re{ det[ p L } reg ) = cos( V ) l[\ n . (7.2) 

n 

Therefore its main advantage is that real rather than complex terms may be summed to 
yield the final result. It still remains a formidable technical problem (although not one of 
principle) to compute the phase rj on the lattice. 
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8 Appendix 



In this appendix we discuss the behavior of Dirac eigenfunctions for large A n and implications 
for the convergence properties of the chiral determinant. 

We begin by considering the two component Dirac eigenvalue equation fl3.15| ), recopied 
here for convenience. 

(<?4 + X n )u n (q) = J d 4 k (iA A (q - k)u n (k) + Ai(q - k)(JiV n {k)) 

(<?4 - X n )v n (q) = / d A k (iA 4 (q - k)v n (k) + A^q - k)(JiU n (k)) . (8.1) 



Consider the limit that X n becomes arbitrarily large (in particular, we want n » Na)- We 
can see immediately that either v n (q) or u n (q) must approach zero in this limit, depending 
on the sign of q^. Let us choose q^ < so that u n (q) is nonzero. The equations in ( |8.1| ) then 
reduce to 

(q 4 + X n )u n (q) = f d 4 k iA A (q - k)u n (k) (8.2) 



up to corrections which vanish as A n — > oo. Note that we can set A^ to zero by a suitable 
gauge transform (which may be different for different solutions due to the special Lorentz 
frame we have chosen above). Thus the eigenfunctions in this limit reduce to free solutions, 
up to a gauge transform. We therefore have 



for arbitrary backgrounds A^ix) which are not necessarily pure gauge. Since the eigenvalues 
of ( |S.1|) are gauge invariant, we also learn that 

Jim A^ = XI (8.4) 



There are several consequences of the results ( |8.3| ) and ( |3.4|) . Consider the Jacobian 
factor 

det{ C L } ■ det[ C R } = det{ C L {C R ) T }. (8.5) 

( p.3| ) implies that for < m,n < N — Na the matrix {C L {C R ) T } mn is close to the 
identity 5 mn , and hence does not contribute to the phase rj. From ( |3.27| ) we know that the 
diagonal corners of the matrix iV — Na < m,n < N are Hermitian. Thus we expect that 
the non-anomalous phase 77 should roughly depend only on the Na x Na submatrices of 
ql,r anc j faofc converges to a well-defined limit as N — > 00 for fixed a. This is intuitively 
plausible, since it is mainly modes with n ^ Na which are affected by the presence of the 
pre-regulated background field. 

The result ( |S.4|) is relevant to the existence and convergence of the limit 

N \A 

lim TT (8.6) 

n /x n 
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where A and A' are different background fields. The existence of ( |8.6| ) requires that very large 
eigenvalues are essentially unperturbed by a sufficiently well-behaved background A^(x) or 
A^x)' . Without this property, the weighting factor of gauge configurations in ( |1 . 1| ) would be 
extremely difficult to compute and could exhibit drastic oscillations due to small changes in 
background field A^x). This question is actually also relevant to the the fermion determinant 
in vector-like models like QCD, and is not specific to chiral models. Unquenched lattice QCD 
computations assume that the ratio flS.6|) converges to its limit for eiqenvalues of order the 
lattice spacing. 

Closer examination of (|8.1|) reveals that large eigenvalues behave as 

^ ~ A° + O(A 2 /X n ). (8.7) 
Using A° ~ (n 1 / 4 /^), where L is the size of our box, we have that 




It is plausible that oscillations in the signs of terms in the sum ( |3.8|) allow it to converge. 
This would imply that ( fj.6[ ) is well-defined. 
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